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The feasibility of using data from the NASA STEREO mission for variable star and asteroseismology studies has been examined. 
A data analysis pipeline has been developed that is able to apply selected algorithms to the entire database of nearly a million 
stars to search for signs of variability. An analysis limited to stars of magnitude 10.5 has been carried out, which has resulted 
in the extraction of 263 eclipsing binaries (EBs), of which 122 are not recorded as such in the Simbad online database. The 
characteristics of the STEREO observations are shown to be extremely well-suited to variable star studies with the ability to 
provide continuous phase coverage for extended periods as well as repeated visits that allow both short and long term variability 
to be observed. This will greatly inform studies of particular stars, such as the pre-cataclysmic variable V471 Tau, as well as entire 
classes of stars, including many forms of rotational variability. The high-precision photometry has also revealed a potentially 
substellar companion to a bright (R = 7.5 mag) nearby star (HD 213597), detected with 5 o significance. This would provide 
a significant contribution to exoplanet research if follow-up observations ascertain the mass to be within the planetary 
domain. Some particularly unusual EBs from the recovered sample are discussed, including a possible reclassification of a 
well-known star as an EB rather than a rotational variable (HR 7355) and several particularly eccentric systems, including 
very long-period EBs. 

Keywords: methods: data analysis - binaries: ecUpsing - space vehicles: STEREO - stars: individual: HD 213597 - stars: 
individual: NSV 7359 - stars: mdividual: V 471 Tau 
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1 INTRODUCTION 

The STEREO mission consists of two satellites in heliocentric 
orbit, each observing the solar environment from the Sun's at- 
mosphere out to the Earth's orbit and together producing 3D 
images of this region (Kaiser et al. 20081. The HI-1 imager on 
the Ahead spacecraft (HI-1 A) has been used to conduct a sur- 
vey for variability of stars down to magnitude 12 in its field of 
view (within about 20 degrees of the Ecliptic Plane, see Figure 
[TJ. The lightcurves produced have been analysed to look for 
exoplanet transits and eclipsing binaries. The value of eclips- 
ing variables in studies of stellar properties and evolution is 
well-known jde Mink et al.|[2(310[ l and the trends from this 
sample, as well as some particular stars of interest, will prove 
informative. The detection of many new eclipsing binaries 
around comparatively bright stars, previously undetected 
in many cases due to shallow eclipses, allows for the de- 
tailed follow-up required to improve upon existing models 
of stellar evolution and behaviour. Also of note is a compar- 
ison of the effectiveness of different algorithms in extracting 
the signatures of eclipsing binaries of different types. 

The rest of this paper is divided into sections giving an 
overview of the STEREOI\i\-\A characteristics, describing 
the data analysis pipeline, the results obtained, prospects for 
further study, the conclusions that can be drawn and the table 
showing all 263 EBs is given at the end. 

Further research is being undertaken to identify and char- 
acterise other kinds of variability, which will be the subject of 
future papers. Some HI-IB data has been extracted (for a re- 
gion between 50 and 60 degrees right ascension and 10 and 
20 degrees declination, which includes a variety of interest- 
ing objects, e.g. V 471 Tau - see Section |4| but this has not 
yet been analysed on a large scale. The HI-IB imager is on 
the side of the spacecraft facing the direction of travel and 
the lightcurves show artefacts resulting from micrometeorite 
impacts, complicating the task of signal extraction. 



with an image bin size of 70 arcsec. A filter limits the spec- 
tral bandpass to 630-730 nm with some contribution around 
400 nm and 950 nm. The characteristics of the CCD on the 
HI-IB are similar, with only slight differences in the spec- 
tral bandpass ( |Eyles et al.||2009^ . Each image is the result 
of 30 summed exposures, each of 40 seconds, with one im- 
age produced every 40 minutes. This prevents saturation by 
cosmic rays, although planetary and cometary incursions do 
cause saturation, and also allows for very faint structures in 
Coronal Mass Ejections (CMEs) and the solar environment to 
be detected. Stars down to 1 2th magnitude are detectable and 
about 650000 stars have been recorded with listed magnitudes 
brighter than 1 1.5, with almost 75000 of these being brighter 
than magnitude 9.5 (Figure [1}. Only the very brightest stars, 
those brighter than about 3rd magnitude, saturate the detec- 
tor but useful lightcurves can nevertheless still be extracted. 
Since these stars are mostly within 20 degrees of the ecliptic 
plane, STEREO/HI- 1 A and HI-IB are therefore searching a 
region of the sky not well observed by dedicated exoplanet 
searches, such as CoRoT jAuvergne, M . et al.'2009 1 {V mag- 
nitudes from 11 to 16), Kepler jKoch et al. 2010 ) (R magni- 
tudes from 7 to 17) and SuperWASP ( |Pollacco et al.|2006> (V 
magnitudes from 9 to 13), and it is also able to view much 
brighter stars than are normally observed. 

Before the data from STEREOmi-lK and HI-IB can be 
analysed, it must be passed through a pipeline (shown in the 
Appendix, Figure [b7) in order to 



• Make a shutterless correction ( Eyles et al.|2009 1 

• Apply a large-scale flat field ( [Bewsher et al.|20l0l l 

• Apply the pointing calibration to determine the positions 
of stars jBrown et al.|2009> 

• Flag out saturated pixels/columns and missing data 
blocks ( Eyles et al. 2009) 

• Carry out a daily background subtraction 



2 STEREOim-\X INSTRUMENTATION AND DATA 
PREPARATION 

The HeUospheric Imager on the STEREO -Pihs?Ld spacecraft 
(HI-I A) provided the data for the large-scale analysis. A sum- 
mary of the data analysis pipeline is given in the Appendix. 
This instrument is described in detail in [Eyles et al.| ( (2009^ 
but for convenience a summary is given here. The HI-1 A has 
a field of view of 20 degrees by 20 degrees and is centred on 
a point 14 degrees away from the centre of the Sun, where the 
solar F-corona is the dominant large-scale structure. The aper- 
ture of the HI-1 A lens system is 16 mm, with a focal length 
of 78 mm. The CCD has a pixel size such that the plate scale 
becomes 35 arcsec on a 2048x2048 chip, with images binned 
2x2 on-board, resulting in 1024x1024 images being returned 



2.1 Shutterless correction 

The STEREOIVH-iA and HI-IB cameras have no shutter and 
remain exposed during the readout process and also during the 
clear process prior to each exposure. Consequently, the top of 
the image is exposed for longer than the bottom of the image, 
resulting in image smearing during readout. The effect of this 
can be expressed as a matrix relationship and a correction for 
the effects of the shutterless operation can be made by a ma- 
trix multiplication. If saturation occurs anywhere in a column, 
then the signal in the affected pixels cannot be corrected in 
this way and the entire column is flagged as bad data. Further 
details can be found in |Eyles et al1 ( |2009l >. 
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2.2 Flat fields 

The sensitivity of the STEREO lYil-l A and HI- IB cameras 
varies spatially in three distinct ways. Firstly, the optical prop- 
erties mean that the solid angle corresponding to each pixel 
may vary across the field of view, leading to large variations 
in scale across the CCD. Secondly, the sensitivity of the CCD 
varies for individual pixels, independently of the sensitivity of 
adjacent pixels. Thirdly, debris and manufacturing variations 
and defects on the CCD generate variations in sensitivity over 
the scale of tens to hundreds of pixels. 

It is unknown to what extent the sensitivity of the CCD 
will vary over time but it is anticipated that the largely sta- 
tionary F-corona in the field of view may be enough to cause 
some degradation. The flat field calibration is updated on a 
regular basis, however, which should prevent this becoming a 
problem. 

The flat field calibrations are performed on STEREOfUl- 
lA and HI-IB images after the removal of the CCD DC bias 
fEyles et al.,2009^ and the application of the shutterless cor- 
rection. 

Images with significant optical ghosts due to saturated 
objects in the data are not used for flat fielding. The method- 
ology is detailed in full in |Bewsher et al.|{2010| l. Observations 
of stars near a region of the frame close to the edge of the solar 
disk are discarded due to additional noise of solar origin. 

2.3 Pointing calibration 

Calibrated instrument pointing solutions have been derived by 
comparing the locations of stars identified in the STEREOfHI- 
lA and HI-IB images with known star positions given in the 
NOMAD astrometric dataset ( |Zacharias et al.||2004l l. In the 
raw (L0.5) data, nominal pointing information is provided, 
however once the data has been processed to LI (by running 
secchi_prep) the calibrated pointing information is provided. 
Full details on the calibration and pointing procedure can be 
found in |Brown etar]p009l >. 

For the STEREO/Hl-IA and HI-IB cameras, the average 
root mean square distance between stars as seen in the image 
and the location given in the catalogue is about 0.1 pixels, i.e. 
7 arcsec. Thus any pointing offsets have a negligible effect on 
the efficacy of the stellar lightcurves. 

2.4 Photometric calibration 

The methodology used to determine the photometric calibra- 
tion is described in full in |Bewsher et al.| 120101, which in- 
cludes more up-to-date flat fields and calibration constants 
than were available for the analysis presented in this paper 
Aperture photometry of well-isolated stars down to 12th mag- 
nitude is used, that have no other stars with a difference of < 1 
magnitude within 1000 arcsec. 




Figure 1. The locations of .stans in the STEREOfHl-lA with mea- 
sured mean magnitude of 9th magnitude and brighter. 

Whilst future large-scale analyses will use the most up- 
to-date calibration of |Bewsher et al.| ( |2010[ >, the photometric 
calibration of the STEREO/Kl-IA used in the analysis pre- 
sented herein is described by the following formula 

, counts , 

STEREO magnitude = dlog( ) (1) 

c 

where d = -2.49501 and c = 18337.7 for STEREO/Hl-IA for 
a 3.5 pixel radius aperture. The number of counts is expressed 
as a flux intensity in units in data numbers, DNi^ ' as provided 
by secchi_prep. For analysing HI-IB, the following formula is 
used (explained more fuUy in |Bewsher etaLlpOTOl l) 

. , ,a* counts, 

STEREO magnitude = d logi ) (2) 

b*c 

where d = —2.5, a = 4 is a correction factor , b = 0.98 and 
c = 97026 are the calibration constants for STEREOfiil-lB 
for a 3.5 pixel radius aperture. 



3 DATA ANALYSIS PIPELINE AND METHODS 

Three methods were used to quantify the variability of each 
star in the STEREO/Hl-lA database with a magnitude of 
10.5 or greater, thus eliminating the faintest stars with the 
lowest signal-to-noise ratio (Figure [T| shows the locations 
of stars with measured mean magnitude of 9th magnitude 
and brighter with more than twenty photometric points on 
their lightcurve). They were: the Box-Least-Squares (BLS) 
I Kovacs et al.|2002 ), the Lomb-Scargle periodogram ( Scargle] 
|1982^ and Stetson's variability index ( Welch & St etson 1993) . 
DDL code to produce a Lomb-Scargle periodogram was ob- 
tained from Armagh Observatory, whilst the others were cus- 
tom implementations. Figure|2]shows the steps in the analysis 
as produced by the programs doing the analysis. 
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Before being passed to the algorithms for analysis, a 7th 
order polynomial fit was carried out for each lightcurve to re- 
move both natural and artificial trends that might obscure the 
faintest transit signals (the top two graphs of Figure [2] show 
an example lightcurve before and after this polynomial de- 
trending). It should be emphasized that whether or not there is 
genuine long term variability that would be removed by this is 
not a concern as the focus is on recovering eclipsing variables 
and exoplanet transits whilst retaining the capacity to detect 
other short term variables, such as RR Lyrae and 5 Scuti stars. 
The IDL built-in routine poly_f it was the basis of this. Each 
lightcurve was then passed through a combined low and high 
pass filter in an attempt to remove some remaining noise. This 
was done by transforming into the frequency domain using the 
IDL built-in routine f f t and applying a mask that blocked all 
features within one standard deviation of the mean before us- 
ing the inverse transform to return to the time domain. This 
process was repeated for each epoch of observations for each 
star and each epoch was analysed separately by the differ- 
ent algorithms (The third graph in Figure [2] shows the BLS 
spectrum). The strongest signals found from analysing all the 
epochs for each star were recorded and stored in an output file 
associated with each field of view, as were a variety of other 
statistics (q.v. Section[B|l. 

The reasons for processing each epoch separately are the 
presence of noise and artefact features that exist within one 
or more epochs but seldom all, such as a planetary incursion 
into the field of view by Venus or Mercury, the existence of 
some features near the beginning or end of some epochs as 
a result of streamers or other solar activity and the presence 
of occasional de-pointing events, typically from micromete- 
orite impacts, all of which would produce a false signal if all 
epochs were analysed together. As the highest values from 
all epochs are stored, these features still give some false sig- 
nals but by analysing data from unaffected epochs a genuine 
value may sometimes supersede the false one. Nevertheless, 
the severity of some noise and artefacts precludes detection of 
faint signals, which is only partially alleviated in some cases 
by a direct visual inspection of lightcurves for signs of gen- 
uine variability. 

Once a star had been identified as an eclipsing binary 
candidate on the basis of the statistics generated, its detrended, 
filtered lightcurve and phase-folded lightcurves on the periods 
of the strongest signals found by the BLS and Lomb-Scargle 
analysis were examined visually to determine if the variability 
detected was likely to be genuine. An initial classification of 
the type of variability was made at this stage. Those that were 
considered likely to be eclipsing binaries were imported into 
the software package Peranso ( jVanmunster|2008} for more 
detailed analysis. There is the possibility that many contact bi- 
nary systems might have been excluded at this stage due to an 
inability to distinguish between different types of variability. 
Indeed, there is some evidence that a significant number of 
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Figure 2. The main stages of the detrending and analysis, demon- 
strated by the star HR 1750. Starting with the top graph showing 
the raw light curve for one epoch of observations, the next graph 
down shows the results of detrending with a 7th order polynomial, the 
third grapli shows the periodogram derived from the detrended 
lightcurve and the bottom graph is the Hghtcurve phase-folded on 
the period of the strongest signal found from the periodogram. 
Note that it is unknown whether the trend removed by the polyno- 
mial fitting is artificial or genuine variability - this potential long term 
variability is not the type of variability being searched for. 

fainter contact binaries may have been excluded in this man- 
ner (Section[6|. 

Within Peranso, a detailed analysis of the BLS spec- 
trum would be used to establish as accurate a period as pos- 
sible using as many epochs as possible, using the software's 
ability to manually exclude observation points that are likely 
to be noise or artefacts. The lightcurve phase-folded from this 
period would be used to finally determine if the star was likely 
to be an eclipsing binary. If it was still considered an eclipsing 
binary, the phase-folded lightcurve would provide estimates 
of the primary eclipse depth and minimum eccentricity. The 
minimum eccentricity would be estimated by visually deter- 
mining the offset of secondary eclipses, where these could be 
distinguished as such, from a phase of 0.5 away from the pri- 
mary eclipse. The estimated errors in the estimate were taken 
as the extremes of the possible ranges of where the secondary 
minimum might be, thus the errors are larger for smoother, 
shallower eclipses as in many contact binaries. The values ob- 
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DISTRIBUTION OF MAGNITUDES OF ECLIPSING BINARIES 



II 



WEIGHTED MEAN MAGNITUDE 



DISTRIBUTION OF PERIODS OF ECLIPSING BINARIES 



PERIOD (DAYS) 



Figure 3. The distribution of the 263 eclipsing binaries observed by 
STEREO lai-l A. by weighted mean magnitude. The fall-off after mag- 
nitude 9 is the result of the decreasing signal quality with magnitude 
increasingly preventing the detection of eclipse signatures. 

tained were then, along with the errors, multiplied by two, 
giving a value that is actually \e x cosco| or effectively a mini- 
mum eccentricity that could be significantly higher depending 
on the value of CO, the argument of periastron. 

The values of the period, primary eclipse depth and |e x 
coscol found from this detailed examination were recorded, 
along with the output of the different algorithms, the star's 
Celestial co-ordinates, its identity (from SiMBAD) and spec- 
tral type (from SiMBAD) as well as the identity of any nearby 
star (within 6 arc-minutes, from SiMBAD) that could poten- 
tially have contaminated the signal and whether or not either 
the star or this neighbour had previously been identified as an 
eclipsing binary. This dataset of 263 eclipsing binaries (Table 
|B1^ was then explored for interesting trends and patterns. 

3.1 Trends 

The distributions of the recorded characteristics show some 
trends that apply to the sample as a whole. Most significantly, 
the distribution of eclipsing binaries by weighted mean mag- 
nitude (Figure [3]( shows a peak between magnitude 8.5 and 
9.0 after which the numbers extracted fall off. This implies 
that, due to a combination of sensitivity and noise, eclipses 
are not being detected around fainter stars since the distribu- 
tion of eclipsing binaries is expected to match the distribution 
of stars in the database, the numbers of which continue to in- 
crease with increasing magnitude. 

The periods given in the statistics are the primary-to- 



Figure 4. The distribution of the 263 eclipsing binaries observed by 
STEREOm-\A by period. 



primary period, although difficulty in distinguishing primary 
from secondary eclipses may cause the resulting value to be 
out by a factor of 2. The distribution of the periods of the 
eclipsing binaries extracted reflects, as expected, the bias due 
to the length of individual epochs (Figure|4]l. In order to get a 
good signal, three transits must be observed within an epoch 
and thus there is a tendency to find periods of less than this 
threshold. Through the presence of secondary eclipses longer 
periods can be found but visual examination is required for 
detecting the very longest period eclipsing binaries, espe- 
cially if the system is eccentric. The star with the longest 
period confirmed of the 263 is HD 72208 with 22.013 days 
(see Section |4] for details), which cannot be found through 
a separate analysis of each epoch but only through all three 
epochs presently available, owing to the significant eccentric- 
ity (|e X coscol = 0.391 ±0.006) and the fact that no more than 
two transits are observed in any epoch. There is a suspected 
eclipsing binary with a potentially longer period showing only 
one primary eclipse in each epoch (HD 173770, see Section 
|4j but a definite period determination cannot be made and it 
was not included in the sample. 

The distribution of the depths of primary eclipses reveals 
a trend that is interpreted as the presence of different popula- 
tions of eclipsing binary (Figure|5j. Only detached and semi- 
detached EBs appear to have a primary eclipse depth greater 
than 0.5 magnitudes, whereas contact binaries predominate at 
depths of 0.3 magnitudes and below. 

The estimates of |e x cosco| reveal that 24 of the 263 
eclipsing binaries have measurable eccentricity (see the Ap- 
pendix for their phase-folded lightcurves). Other EBs in the 
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DISTRIBUTION OF PRIMARY ECLIPSE DEPTHS 



DISTRIBUTION OF ECCENTRICITIES WITH ERROR MEASUREMENT 
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Figure 5. The distribution of the 263 eclipsing binaries observed by 
STEREOIHl-lA by primary eclipse depth. 



sample may have eccentricity, as the value of the argument of 
periastron, to, is unknown and there are some that show a hint 
of eccentricity but the margin of error does not preclude the 
possibility that they might have zero eccentricity, including 
a couple that are exactly on the line where the error margin 
equals the estimate of |e x cosa)|. Thus a small but significant 
fraction of the sample are eccentric, with more than 5% show- 
ing sigificant eccentricity. The sample of eccentric systems is 
not quite large enough to show statistically valid trends, unfor- 
tunately, although the indications are that small eccentricities 
dominate with larger eccentricities being widely distributed. 
Nevertheless, this result is broadly in agreement with recent 
findings from the Kepler mission ( |Prsa et al.|2010| l, however 
the sample of stars observed by STEREOfHI-lA is not limited 
to main sequence stars or by other selection criteria intended 
to maximise the chances of finding an exoplanet transit. It is 
therefore a useful result that it appears to find a similar distri- 
bution of |e X cosco| to the Kepler sample as it implies that 
eccentricity is not dependent upon spectral type. Eccentric 
systems are expected to circularise over time and thus these 
systems are much-studied in the field of stellar evolution to 
model interactions between multiple stars as an interaction 
with a third body is usually required for them to be formed 
jClaret & Willems|2002| l. 

The distribution of spectral types of the eclipsing binary 
primary stars shows a predominance of younger, hotter stars 
of earlier spectral types (Figure[8j. The main limitation of this 
distribution is that the spectral type of the closest star to the 
co-ordinates of the perceived eclipsing binary may not always 
be the star that is actually eclipsing, due to the low spatial 



Figure 6. The estimates of the minimum eccentricities (|e x cosQ)|) of 
the 263 eclipsing binaries observed by STEREO/Rl-IA, where mea- 
surable, against the errors for each measurement. Each open circle 
represents a single eclipsing binary. Those points above the line rep- 
resent eclipsing binaries for which an eccentricity of zero is outside 
the range of the errors in the measurement. 



DISTRIBUTION OF ECCENTRICITIES AS |e*cos(w)| 



Figure 7. The estimates of the minimum eccentricities (|e x cos(a|) of 
the 263 eclipsing binaries observed by STEREOfHl-lA, where mea- 
surable. 
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RELATIONSHIP BETWEEN PERIOD AND SPECTRAL TYPE 



DISTRIBUTION OF BLS SIGNAL DETECTION EFFICIENCY 
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Figure 8. The 232 eclipsing binaries observed by STEREOIli\-\K 
witli a known spectral type, showing the relationship between spec- 
tral type and period. The most common spectral type is A2, with 22 
stars in the sample. There ai'e no clear patterns connecting the spec- 
tral type with period, however there are some cases where the specific 
star causing the eclipses has not been identified and the spectral type 
is of the star closest to the target co-ordinates (see AppendixjAj. 



resolution and the possible contamination of lightcurves by 
nearby stars. The distribution will also be affected by any bias 
in the STEREOIYH-IA. database, with only stars brighter than 
magnitude 10.5 being analysed en masse (although one or two 
fainter examples were identified visually). 

The distribution of eclipsing binaries by their signal de- 
tection efficiency (SDE) from the BLS algorithm (Figure |9j 
shows two indistinct peaks, interpreted as corresponding to 
two different types of eclipsing binary. The peak with the 
highest SDE corresponds mostly to detached binaries, whose 
lightcurves are a closer match to the box-like shape the BLS 
algorithm checks for. The broader peak with the lower SDE 
corresponds to contact binaries, which have smoother, more 
sinusoidal lightcurves. An Algol type binary might be the 
highest-rated by the BLS algorithm in its field but contact bi- 
naries can sometimes be near the bottom of such a list and thus 
are not reliably detected by this algorithm. The Lomb-Scargle 
periodogram, (distribution shown in Figure[TO]l, preferentially 
reveals contact binaries, although it typically does give a sig- 
nal for many detached binaries as well, albeit increasingly less 



BLS SIGNAL DETECTION EFFICIENCY 

Figure 9. The distribution of the 263 eclipsing binaries observed by 
STEREO/Hl-IA by BLS signal detection efficiency. 



so for shallower, more box-like eclipses (Figure [T7). Com- 
bining the BLS and the Lomb-Scargle methods thus produces 
a positive signal for the vast majority of the eclipsing bina- 
ries, with those rated poorly by both being stars with fewer 
data points (both algorithms give higher scores for more data 
points, thus stars with less data available are rated less highly), 
those with very long periods and those with comparatively low 
signal-to-noise that were detected by visual examination. 

The effect of the period on the output of the BLS and 
Lomb-Scargle methods was also investigated (Figures[T2|and 
\T3\ . The dependence discovered is again likely representative 
of the shape of the transits. 

As a variety of noise signatures were encountered in the 
dataset, some very clearly obscuring genuine signals, produc- 
ing both false positives and false negatives, an attempt to eval- 
uate the effects of noise was made (Figures [l4][T5][T6]and[l7|. 
Together these appear to show a tendency to overlook shallow 
eclipses. This may be partly due to a mixture of red and white 
noise and also due to problems classifying contact binaries as 
such when their lightcurves are not clear enough to distin- 
guish them from sinusoidal variations. Blending also has 
a tendency to obscure faint signals, with the brightest ob- 
ject in the photometric aperture dominating but showing 
additional signals overlaid, or with stars of similar mag- 
nitudes showing constructive or destructive interference 
resembling an amplitude modulation. This has the effect 
of lowering, broadening and sometimes splitting peaks in 
a periodogram analysis, which affects the ability of an al- 
gorithm to identify a low-amplitude signal. 

Since a large proportion of the sample of eclipsing bi- 



© 0000 RAS, MNRAS 000, 000-000 



8 K.T. Wraighty Glenn J. White, D. Bewsher and A. J. Norton 



DISTRIBUTION OF PEAK LOMB-SCARGLE POWER 
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Figure 10. The distribution of tlie 263 eclipsing binaries observed by 
STEREOIHl-lA by peak power from a Lomb-Scargle periodogram 
analysis. 



Figure 12. The dependence of the BLS algorithm on the period of a 
signal. The influence is believed to be the shape of the light curves of 
variables with increasingly short periods, rather than the period itself. 
Each open circle represents a single eclipsing binary. 



DIFFERENT TYPES OF ECLIPSING BINARIES 
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PEAK LOMB-SCARGLE POWER AS A FUNCTION OF PERIOD 
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Figure 11. Comparison of the ability of the BLS and Lomb-Scargle 
algorithms to extract the signals of eclipsing binaries. Each eclipsing 
binary is represented by a circle with a radius proportional to the depth 
of the primary eclipse. Points above the horizontal line have a SDE 
high enough to be considered a detection by the BLS algorithm, whilst 
those to the right of the vertical line may similarly be considered to 
have been detected by the Lomb-Scargle algorithm, in the absence of 
artefacts and noise producing large numbers of similar signals. Those 
rated poorly by both were typically found through visual examination. 
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Figure 13. The dependence of the Lomb-Scargle method on the pe- 
riod of a signal. Each open circle represents a single eclipsing binary. 



naries appear to be previously undetected or, at least, uncon- 
firmed as such by previous observations, a search was made 
for characteristics in common to understand what, if anything, 
might be different about the new cases that might explain why 
they had previously gone undetected (Figures [Tsji. As might 
be expected, there is a tendency for the new cases to have shal- 
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Figure 14. The dependence of primary eclipse depth on the Stetson 
index. The sharp cut-off at about -0.5 shows that eclipses with a pri- 
mary eclipse depth less than 0.3 magnitudes are considerably more 
difficult to detect through noise or variability. Each open circle repre- 
sents a single eclipsing binary. 



NOISE BOUNDARY IN DETECTING ECLIPSING BINARIES 




LOG10 STETSON INDEX 

Figure 15. The dependence of peak power from a Lomb-Scargle pe- 
riodogram analysis as a function of Stetson index. The cut-off point at 
about -0.5 shows that variability or noise above a certain level is very 
effective at obscuring the signals of eclipsing binaries, regardless of 
their peak power. Each open circle represents a single eclipsing bi- 
nary. 



Figure 16. The dependence of signal detection efficiency as a func- 
tion of Stetson index. The cut-off point at about -0.5 shows that vari- 
ability or noise above a certain level is obscuring the signals of some 
eclipsing binaries, although the apparent clustering of signals around 
a moderate value and especially the higher ratings for deeper eclipses 
shows that the Stetson index is also to some extent reflecting the vari- 
ability due to eclipses. Each eclipsing binary is represented by a circle 
with a radius proportional to the depth of the primary eclipse. 



PRIMARY ECLIPSE DEPTH AS A FUNCTION OF MAGNITUDE 
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Figure 17. The dependence of primary eclipse depth as a function 
of weighted mean magnitude. After about magnitude 8.5 shallower 
eclipses are not detected, with the threshold rising nearly linearly. 
Each open circle represents a single eclipsing binary. 
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COMPARISON OF KNOWN AND NEW ECLIPSING BINARIES 
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Figure 18. Comparison of the distribution of known and new eclips- 
ing binaries with respect to the depth of their primary eclipses. The 
tendency of the newly-discovered systems to show fainter eclipses re- 
flects both the sensitivity of the STEREO fRl-l A and the expected bias 
of previous surveys in discovering deeper eclipses. The centre box 
contains those EBs in the second and third quaililes with the median 
represented by a thick line within, whilst the outer bars represent an 
approximately 95% confidence level and outliers are individual open 
circles. 



lower eclipses than those previously known. There does not 
appear to be any significant difference in either the periods or 
magnitudes found, however. 

A previously known tendency for eclipsing binaries with 
long periods to be more likely to be eccentric iMat hieu &| 
|Mazeh 1988) is confirmed by the STEREOimAA observa- 
tions (Figure[T9]l. 

To summarise briefly, the BLS is the algorithm of choice 
for extracting box-like eclipses, e.g. Algol binaries or exo- 
planet transits. The Lomb-Scargle algorithm is best at extract- 
ing sinusoidal-like signals but also has some limited capabil- 
ity to extract other regular signals. The Stetson index is able 
to extract large amplitude variables of all types. 



4 RESULTS 

Of the 263 eclipsing binaries extracted, 122 are not recorded 
as such in SiMBAD. Some of these are, however, known vari- 
ables of other types or are suspected of variability. Owing 
to the purely photometric data available from STEREOfHI- 
lA, there remains the possibility of mis-classification, how- 
ever care was taken to include only stars whose phase-folded 
lightcurves appeared to be clear cases. Thanks to the excel- 



Figure 19. The dependence of the minimum eccentricity of the 263 
eclipsing binaries observed by STEREO/Hl-IA, as measured by \e x 
cos(a|, on the period of eclipses. Each open circle represents a single 
eclipsing binary. 



lent phase coverage resulting from the STEREOfHI-lA ob- 
servational cadence of 40 minutes maintained for up to 19.44 
days, with three or four such epochs for most stars, the new 
observations are in many cases clearer than those previously 
available and in some cases the existing classification may be 
in need of revision. A total of 24 of the eclipsing binaries have 
measurable eccentricity and their lightcurves are shown in the 
Appendix (Figure [52^. 



4.1 Individual stars of interest 

4.1.1 HD 213597: substellar transit candidate 

The best potentially substellar transiting candidate is HD 
213597. This star is recorded in SiMBAD as being of spectral 
type FO and is not suspected of any variability. The observed 
eclipses are box-like, with very sudden ingress and egress 
phases, and are uniformly about 25 mmag deep, with a pe- 
riod of 2.4238 days between successive eclipses (Figure [20)1. 
Secondary eclipses are not evident. 



4.1.2 V 471 Tau: cataclysmic variable progenitor 

Although it does not feature in the sample, as the eclipses 
by the white dwarf secondary are too short to be observed, 
V471 Tau nevertheless is an important star and the STEREO 
observations in both HI-IA (Figure [22]l and HI- IB (Figure 
[23| > have the potential to inform studies of the magnetic activ- 
ity cycle of the red dwarf primary (eg. ( [O'Brien et al.|200T] ), 
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Figure 20. Lightcurve of the star HD 213597, phase-folded on a pe- 
riod of 4.8476 days. 

jHussain et al.|2006l l and ( [Kaminski et al.|2007^ ). As a cat- 
aclysmic variable progenitor, with an unexpectedly massive 
red dwarf primary and the most massive white dwarf in the 
Hyades ( [O'Brien et al.|2001| ), studies of this star may have a 
considerable effect on understanding the mechanism behind 
type la supernovae, thereby influencing the calibration of the 
extragalactic distance scale and much of modem cosmology. 
Note that the spectral bandpass of the Hl-IA and HI-IB im- 
agers both allow a component of blue light through ( [Bewsher] 
|et al.|201^ and the white dwarf secondary is estimated to con- 
tribute 6 % of the flux in this part of the spectrum ('O'Brien' 
|et al.| |2001|l. Similar behaviour can be seen in STEREO ob- 
servations of RS CVn variables, e.g. SZ Psc, and shows that 
STEREO is a useful resource for studies of similar variability 
that require photometry with both long and short baselines. 

4.1.3 Highly eccentric eclipsing binaries 

The most eccentric undiscovered eclipsing binary in the sam- 
ple is HIP 92307, with \e x cosQ)| = 0.500 ±0.020 (Figure 
\24) . Both primary and secondary eclipses are about 30 mmag 
deep, indeed it is not clear which are the primary and sec- 
ondary eclipses. This star is listed as spectral type A2V. 

The longest period eclipsing binary in the sample is 
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Figure 21. Close-up of the transit of HD 213597, phase-folded on 
a period of 2.4238 days. This lightcurve was constructed using the 
latest data available, including updated flat-fields as well as data 
from STEREO/Ul-IB. 



HD 72208, with 22.0130 days between primary eclipses and 
\e X coscol = 0.391 ± 0.006 , making it one of the more ec- 
centric as well (Figure [25). This star, spectral type B9p, is a 
spectroscopic binary. Both primary and secondary eclipses are 
about 60 mmag deep. 

Possibly the longest period eclipsing binary so far ob- 
served in the STEREO/Rl-IA data may be HD 173770. As 
there are multiple possible periods that can fit the lightcurve, it 
was not included in the sample. The best-fitting period that has 
been found with the data available shows a period of 25.145 
days, which implies a minimum eccentricity \e x cosco| = 
0.108 ±0.012. The star shows some intrinsic regular variabil- 
ity as well, although the exact type has not been classified. 
The host star is listed in SiMBAD as being of spectral type 
B6V and is a suspected triple system. With the separation given 
for these components, it is unlikely either are responsible for 
the eclipses and this may instead be a quadruple system. The 
lightcurve is included in Figure [B2[ for comparison with the 
other measurably eccentric eclipsing binaries. 
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Figure 22. Lightcurve of the star V471 Tau, as observed by the 
STEREO HI-IA imager 3267 is the catalogue number of the star in 
the STEREO HI-IA field of view for right ascension: 50-60 degrees 
and declination: 10-20 degrees. 

4.1.4 NSV 7359: (3 Cepheid showing the Kozai effect 

A well-observed bright star, NSV 7359 shows relatively deep 
eclipses with a period of 9.1999 days that have not been 
recorded in SiMBAD (Figure |26[ l. This star is a known |3 
Cepheid, with an amplitude listed as 30 mmag. Given the 
volume of observations, it is very surprising this variability 
has not previously been recorded. It is in a fairly crowded 
area with many stars close by in the STEREOI^-\h field of 
view, yet none of these are sufficiently bright to contaminate 
its lightcurve to the extent of eclipses as deep as those ob- 
served. The eclipses are therefore genuine and should help to 
further inform the parameters previously determined for this 
star. Interestingly, this star is also shown in SiMBAD to be 
a spectroscopic binary with a period of 0.2872 days, which 
would mean the eclipsing companion is a tertiary component 
and not the secondary. The Kozai effect ( |Kozai||T962| l may 
have caused the inclination of the system to change so that 
eclipses are now visible. 

4.1.5 HR 7355: rotational variable or contact binary? 

One of the better-observed stars in the sample is HR 7355 
(a.k.a. HIP 95408 and HD 182180). STEREOim-iA shows 



Figure 23. Lightcurve of the star V471 Tau, as observed by the 
STEREO HI- IB imager. 3130 is the catalogue number of the star in 
the STEREO HI-IB field of view for right ascension: 50-60 degrees 
and declination: 10-20 degrees. The scale of the y-axis has been ad- 
justed to tit the data range. 

what appears to be a clear-cut case for an eclipsing binary as 
the phase-folding on 0.5214 days indicates in Figure|27] Other 
studies have noted this periodicity and associated it with a ro- 
tation period, however, rather than attributing it to a secondary 
companion, even the most recent which show the variability 
to be eclipse-like (Oksala et al.||2010l >. We suggest that the 
photometry from STEREO is enough to warrant a re-appraisal 
of the nature of this system. It has been left in the sample as 
an eclipsing binary since from the photometry alone that is 
what it appears to be. In the wider context it is clearly an un- 
usual and fascinating, chemically peculiar star that continues 
to attract further study. 

4.2 False positives 

A small number of stars showing variability on a similar scale 
to that expected of an exoplanet or brown dwarf transit could 
be ruled out as such from photometric data alone, or in com- 
bination with other information known about the host star. 

Total eclipses just 30mmag deep were observed 
around either HIP 248 or HIP 247 (near-identical sig- 
nals and in the same or adjacent pixels in STEREOfHl- 
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Figure 24. Lightcurve of the star HIP 92307, phase-folded on a period 
of 4.7911 days. |excos(B| = 0.500±0.020. 
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Figure 25. Lightcurve of the star HD 72208, phase-folded on a period 
of 22.0130 days. |excosa)| = 0.391 ±0.006. 



lA). These stars are both of spectral type FO, however 
the eclipses include a substantial amount of time in the 
ingress and egress phases. The eclipses are consistently 
about 30 mmag deep and 2.2604 days apart (Figure [21). 
This is a transit by a larger, stellar body rather than a 
brown dwarf but a low-mass star remains a possibility. 

The star HD 222891 has similar features to HIP 248, 
however the eclipses are also too deep, at 40 mmag, with 
the host star of spectral type F8, to be due to a brown 
dwarf companion. It may nevertheless be a low-mass star of 
some interest. 

Of all eclipse-like features extracted from the 
STEREOfUl-lA database, the shallowest primary eclipses 
were 7 mmag deep around the star 18 Sgr. This lightcurve ap- 
pears more like a contact binary, although the raw lightcurve 
is also somewhat irregular and the behaviour might instead 
be due to rotational microvariability. It was not considered 
conclusive enough to be labelled as an eclipsing binary and 
does not feature in the sample, although more data may help 
to clarify the lightcurve and its classification. 

NSV 1321 is the most marginal detection in the sam- 
ple. In order to verify whether the 10 mmag eclipses ob- 
served were genuine, differential photometry was carried 
out, however this appeared to introduce variability from a 
neighbouring star - it is located in a crowded field with several 



stars showing clear variability and others blended. Although 
the star itself is classified as variable, the NSV catalogue as 
shown in SiMBAD indicates this variability was considered 
doubtful by the catalogue's compilers. With the host star of 
spectral type B8V this would not be an exoplanet signal and 
is most likely a grazing eclipse or stellar companion. Al- 
though the possibility of microvariability or a blend with a 
neighbouring star cannot be wholly ruled out, none of the 
neighbouring stars show a periodicity of 2.2663 days so this 
was included in the sample. 

Another star showing shallow eclipses that was ruled 
out from being a brown dwarf is II Cnc (period from 
STEREO[m-\A of 1.0684 days, eclipses 25 mmag deep). 
This is listed in SiMBAD as a BY Draconis type variable of 
spectral type G8V, however there is no mention of it show- 
ing eclipses. The shape of the lightcurve is indicative of ei- 
ther a contact binary or grazing stellar transits with secondary 
eclipses present. The period is more representative of a con- 
tact binary. The secondary may nevertheless be a low-mass 
star, making this a relatively unusual system for a contact bi- 
nary. 
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Figure 26. Lightcurve of the star NSV 7359. phase-folded on a period 
of 9. 1999 days. 
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Figure 27. Lightcurve of the star HR 7355, phase-folded on a period 
of 0.52 14 days. 



5 SURVEY PROSPECTS 

Although the analysis thus far done has been restricted to re- 
covering eclipsing binaries and searching for exoplanets, the 
prospects for recovering variables of different types has been 
investigated. 

The prospects for exoplanet extraction are summarised 
in Figure|29]for the declination strips away from the Galactic 
Plane for the ranges of periods and magnitudes shown. This 
indicates reasonable prospects for finding a signal, if any ex- 
ist around stars of about 7th magnitude or brighter but with 
a sharp fall-off for fainter stars, such that the probability of 
extraction is negligible for stars fainter than 8th magnitude 
from STEREOfHI-lA data alone. Note that the most signifi- 
cant substellar candidate so far extracted is magnitude 7.5. 
For reference, the signal of HIP 247 and HIP 248 (magnitude 
7.5, period 2.2604 days) is 30 mmag deep detected with a sig- 
nificance of 5 a and the 30 mmag deep signal of HR 1750 
(magnitude 6.5, period 3.3150 days) is detected with a signif- 
icance of 8 o. With the inclusion of STEREO/Hl-IB data the 
number of data points will more than double (this satellite is in 
an Earth-trailing orbit and stars remain in the field of view for 
longer) and the significance of extracted signals will increase, 
as will the potential for extracting faint signals. 

For eclipsing binaries, Figures p4l|17| and|28|indicate that 



more may yet be found and Figure[TT|shows that the combina- 
tion of the BLS and Lomb-Scargle methods is a good way of 
extracting their signals. The difficulty here is that, from photo- 
metric data alone, it has not been possible to conclusively clas- 
sify them as eclipsing binaries. This problem largely affects 
contact binaries, which have lightcurves that might easily be 
confused with elliptical variables and short period pulsating 
stars. The ability of the Lomb-Scargle method to extract con- 
tact binaries suggests it would be a useful tool for extracting 
other variables also. 

Asteroseismological studies are likely to require that 
low-level variability be detected. If a signal of 10 mmag, com- 
parable to an exoplanet transit, is required, then Figure [29] 
might also be an indicator of the prospects for extracting mi- 
crovariable behaviour, however this would assume that the al- 
gorithm chosen to extract the microvariable signal is equally 
effective at discerning faint signals as the BLS. Nevertheless, 
the prospects would appear good for the brighter stars (7th 
magnitude or brighter) and these are also the most amenable 
to high-resolution spectroscopic follow-up by other observa- 
tories. Investigations of the signals extracted by the Lomb- 
Scargle method show that the presence of systematic sources 
of noise inhibits the detection of a test signal, however, so 
additional cleaning would be required. Whilst the observ- 
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ing cadence of 40 minutes will impair detailed studies of 
very short period variables, the ongoing examination of 
the data that is including STE/fEO/HI-lB observations is 
nevertheless able to obtain very clear photometry from a 
number of known 5 Scuti stars (e.g. DX Cet) and poten- 
tially detect new variables down to magnitude 11. 



6 CONCLUSIONS 

The STEREOfHI-lA observations are clearly a valuable re- 
source suitable for a wide range of variability studies. With 
263 eclipsing binaries extracted from the database, 122 of 
which have not previously been identified as such, there is 
the potential for the discovery of a number of new and in- 
teresting objects. The area of the sky being covered and the 
brightness of the stars being observed are features unique 
to STEREOfHl-lA and simultaneously allow for high qual- 
ity follow-up observations of interesting objects. There is the 
potential to discover transiting exoplanets (Figure [29]), how- 
ever this is at the very limit of the survey's sensitivity and 
would only be possible for the stars with the cleanest sig- 
nals. A 30 mmag transit around a magnitude 7.5 star is de- 
tectable with 5 a significance and a 30 mmag transit of a 
magnitude 6.5 star is detectable with 8 O significance from 
STEREOfHI-lA data alone. With the inclusion of data from 
the second satellite, STEREO/Hl-IB, this will improve, e.g. 
from STEREO fHl-l A data alone the 50 mmag transits of the 
magnitude 8.5 star BD +03 263p are detected with 7 a sig- 
nificance but with the inclusion of STEREOfHI-lB data they 
are detected with nearly 13 O significance. 

The analysis presented here may have some impact on 
stellar evolution studies, in particular relating to the propor- 
tion and distribution of eccentric eclipsing binaries. As the 
newly-discovered eclipsing binaries are comparatively bright, 
the masses, radii and effective temperatures can be determined 
more accurately in follow-up observations and can therefore 
provide more stringent tests of stellar evolution models (e.g. 
^Torres e t al.||2006[l, ([VandenBerg eraLl [2006), ( |Claret &| 
[Willems 2002|i and jStassun et al.'2006)). The new eclipsing 
binaries will also help to produce more accurate distances to 
their host stars and any stars those are associated with jSouth-| 
[worth et al.|2005| l. 

The potential to improve upon existing phase coverage, 
especially relating to that obtained during eclipses, and times 
of maxima and minima in a variability cycle can help to refine 
parameters for these important stars. In some cases, the capa- 
bility of STEREO to deliver good photometry over timescales 
of days and weeks in a single pass, combined with repeated 
observations after about a year, are ideal for observing long- 
term changes in variability. This is of particular importance in 
studies of magnetic activity cycles (eg. V471 Tau) and partic- 
ular classes of variables, such as RS CVn stars (eg. SZ Psc). 
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Figure 28. The proportion of stars of different magnitudes that have 
been found to be eclipsing. The small number of stars brighter than 
magnitude 4.5 in the STEREO/Hl-IA field of view is responsible 
for the shape at the bright end of this graph, with only three EBs 
observed. 



The eclipsing binaries newly discovered by this analysis 
show a tendency to possess shallower transits than those pre- 
viously known (Figure [Tsj, partly due to the quality and ca- 
dence of the photometry. With many of the host stars being of 
early spectral types, late B and early A, it is possible that many 
of these systems may be comparatively young (Figure [Sf. A 
small but significant proportion are in eccentric orbits (Fig- 
ures[6|and[7] l, supporting recent observations by Kepler jPrsa| 
|et al!pOfo[ rSome particularly interesting systems have been 
identified for potential follow-up observations, including one 
good brown dwarf candidate (HD 213597, Figure [20} . The 
proportion of stars of different magnitudes found to be eclips- 
ing indicates that more may yet be found as Figure [28[ shows 
decreasing numbers observed with fainter stars, whereas the 
real number should be a constant proportion of the popula- 
tion. 

Evaluation of the algorithms used to extract eclipsing bi- 
naries suggests that, if exoplanets alone are the target, then the 
Box-Least-Squares (BLS) algorithm is capable of effectively 
extracting them (Figure [29}. For a more general search for 
eclipsing binaries, however, the combination of the BLS along 
with the Lomb-Scargle periodogram analysis will be required 
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Figure 29. The percentage of stars where a 10 mmag test echpse su- 
perimposed on the real lightcurve was successfully returned as the 
strongest signal by the BLS algorithm. Each epoch of each star was 
examined separately and only those that returned the test signal within 
0.04 days were considered successes. Detection of a harmonic was 
not considered a success. This graph thus represents the chance of 
detecting a transiting 'Hot Jupiter' for stars of a given magnitude for 
the range of periods shown. The stars tested were those with a dec- 
lination between -10 degrees and 40 degrees and having a weighted 
mean magnitude within 0.05 magnitudes of the points shown. The 
trial periods ranged from 0.5 to 9.5 days in half-day steps. The dip 
in success rates at 6.5 days is the result of many stars having parts 
of their lightcurve masked out to avoid contamination by solar activ- 
ity, whilst the dip at 9.5 days is the result of fewer stars remaining 
in the field of view for long enough for a period of this length to be 
checked. The small increases in success rates for fainter stars at these 
periods are due to artificial trends associated with the length of an 
epoch, which also contributes to an increase in the number of false 
positives for brighter stars although this number is small compared 
to the decreases observed for the above reasons. The best substellar 
candidate so far extracted has a period of 2.4238 days and a host 
star of magnitude 7.5: a 25 mmag transit is detectable with 5 o 
significance at this level. 

in order to extract a high proportion of these objects (Figure 
[TTJ. The Stetson index can be used as a measure of noise, as 
well as variabihty, ahhough it does rate some Algol type bi- 
naries very highly so extreme values are not necessarily false 
positives (Figure [14). 

Further research being done includes the development 



of a matched filter algorithm and follow-up observations of 
key targets has begun. A further two eclipsing binaries show- 
ing eccentricity have been discovered (TYC 1422-1328-1, 
magnitude in STEREOimAK of 9.5, period 3.0999 days, 
\e X cos(b| = 0.175 ± 0.025 with primary and secondary 
eclipses both 0.1 mag deep and BD -1-09 485, magnitude in 
STEREOim-\A of 9.15, period 4.5293 days, \e x cosco| = 
0.125±0.025 with primary eclipses 0.15 mag deep and sec- 
ondary eclipses 0.1 mag deep), amongst numerous other 
short-period variables. Periods have been extracted for a 
variety of stars showing sinusoidal-like variations in their 
lightcurve and some stars showing changes in behaviour on 
a year-by-year basis have been observed (eg. V 471 Tau, SZ 
Psc). Additional findings will be reported in due course. 
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APPENDIX A: ECCENTRIC ECLIPSING BINARIES 

The measurably eccentric eclipsing binaries are shown in Fig- 
ure |B2] along with the lightcurve of HD 173770, which al- 
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though not included in the sample owing to an inability to 
uniquely determine its period due to the small number of 
eclipses observed, nevertheless appears to be an eccentric 
eclipsing binary. The \e x cos(X)| measured for this star is 
0. 108 ± 0.012 and the lightcurve shown is a phase-folding on 
a period of 25.145 days. The values of \e x cos(X)| measured 
for the other stars can be found in Table IB 1 1 The identities of 
the stars are as given in the table, however in the cases where 
nearby stars may have contaminated the signal or be the real 
source of eclipses it was often the case that the eclipses would 
feature in multiple lightcurves and the ones selected for anal- 
ysis were those that were clearer and had deeper eclipses and 
the higher SDE. Where this occurred, the table gives the name 
of a star that may be contaminating the lightcurve or be the 
real source of eclipses. 



APPENDIX B: TABLES AND FIGURES 

The sample of 263 eclipsing binaries, with the various statis- 
tics that have been collected for each is given here. Where 
data is not available, for example because secondary eclipses 
were not visible or the spectral type has yet to be determined, 
the entry is given as NA. Right ascension and declination are 
given in degrees. Primary eclipse depth and weighted mean 
magnitude are both in units of magnitudes. Periods are in 
units of days. The Box-Least-Squares (BLS) signal detection 
efficiency (SDE), the peak power from a Lomb-Scargle pe- 
riodogram analysis and the Stetson variability index are all 
arbitrary scalar units. The estimates of minimum eccentric- 
ity (|e X cosco|) and the estimated errors in those estimates 
are all scalar values. The spectral type given is for the closest 
star to the listed co-ordinates from SiMBAD, the identity of 
which is given in the next column. The alternative nearby star 
is the identity of a nearby star that might be the source of the 
eclipses, rather than the closest star (this frequently occurs due 
to the 70 arc-second per pixel resolution of STEREO A, 
as well as being due to nearby bright eclipsing binaries). The 
Known EB is a binary flag set to 1 if either star is recorded as 
showing eclipses and otherwise. 
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Figure Bl. The STEREOI¥A-\K data analysis pipeline. 
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Figure B2. Lightcurves of all measurably eccentric EBs in the sample and the best-fitting lightcurve of HD 173770 (see Section|4j. 
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Figure B3. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B4. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B5. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B6. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B7. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B8. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure B9. Lightcurves of all EBs in the sample including HD 173770 (see Section|4j. For Supporting Information only. 
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Figure BIO. Lightcurves of all EBs in the sample including HD 173770 (see Section[4j. For Supporting Information only. 
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Figure Bll. Lightcurves of all EBs in the sample including HD 173770 (see Section[4|. For Supporting Information only. 
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Figure B12. Lightcurves of all EBs in the sample including HD 173770 (see Section[4|. For Supporting Information only. 
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Figure B13. Lightcurves of all EBs in the sample including HD 173770 (see Section[4|. For Supporting Information only. 
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NA 
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6.5489 
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NA 


NA 


F2 
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NA 


1 
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28.574 
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NA 
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0.06538 


109.62 


NA 
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8.6189 


6.2292 
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52.623 


NA 


NA 


F8 
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NA 


1 
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12.8988 
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1.9631 


7.9653 


5.3368 
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0.015 


B9 
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NA 
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3.7474 
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F8 
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NA 


1 
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NA 


NA 


A2 
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B3V 
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NA 
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NA 


NA 


AO 
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41.01 


0.005 
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V818Tau 
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0.19 
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5.9208 
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NA 


1 
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0.16955 
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13.2113 
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B3V 
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NA 
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13.6861 


0.13 
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NA 


NA 


KO 
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NA 
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FO 
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19.976 
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A3 
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NA 
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14.7035 


0.3 


4.4507 


9.0733 


5.7337 
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2.2604 
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NA 


NA 
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NA 


NA 
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A5 
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10.246 
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32.464 


3.76935 


0.18 


2.3682 


6.6093 


8.0387 


0.61444 


57.542 


0.005 


0.015 


A2 


DPCet 


NA 


1 


46.2098 


8.11422 


0.29 


5.9971 


9.3981 


6.365 


0.050353 


41.737 


0.13 


0.01 


F8 


HIP 14334 


NA 





51.7923 


9.73267 


0.07 


3.5732 


3.7049 


7.3388 


10.227 


38.366 


NA 


NA 


B9Vn 


ksiTau 


NA 


1 


121.911 
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0.06 


22.013 


6.6731 


4.5847 


0.083666 
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0.006 
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NA 





133.069 
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0.4 
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3.5363 
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64.636 


NA 


NA 


A2 


TUCnc 


NA 
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143.758 
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3.0095 


9.009 
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0.005 


0.04 


A2 
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NA 
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0.08 
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8.3795 


3.1369 


0.0061672 


42.305 


NA 


NA 


F8 
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2.445 


5.5845 
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0.03 


A2V 


TXLeo 


NA 
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162.333 
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6.7082 
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32.705 
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NA 
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NA 
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NA 
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NA 
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0.06 
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F2 
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NA 
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NA 


NA 
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NA 
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0.039807 
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NA 


TYC 4913 
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0.3 


0.4434 
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3.3536 
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164.775 

166.176 

169.152 

171.997 

173.589 
181.587 

192.411 

197.118 

208.716 

220.288 

225.243 

323.766 
328.78 

334.768 

338.933 

345.33 
349.842 

356.826 

356.162 

190.304 

196.428 
197.823 



3.54405 

2.97221 

1.90306 

1.92144 
-5.5271 

4.91402 

6.07913 

2.68878 

2.38036 

5.09566 

8.51894 
-3.7349 

9.79926 

2.64167 

0.692436 
-6.4376 

8.87043 

8.08669 

8.84879 



13.0246 
-12.37 

11.1059 



0.3 

0.075 

0.3 

0.18 

0.06 
0.4 

0.5 

0.2 

0.42 

0.18 

0.85 

0.06 
0.25 

0.5 

0.7 

0.6 
0.25 

0.28 

0.04 

0.13 

0.1 
0.22 



1.0713 

1.0847 

0.6288 

1.0143 

2.5670 
2.7113 

1.3086 

2.7323 

6.6950 

0.3846 

2.3273 

2.4449 
1.8636 

2.1597 

0.556 

2.8536 
3.03095 

0.4814 

1.59495 

3.14445 

1.581 
2.99275 



9.5849 
8.591 

9.8227 

8.8186 

6.5105 
10.121 

9.1087 

7.4267 

8.4903 

9.4452 

4.9170 

7.5587 
9.6842 

9.9316 

10.003 

9.4475 
9.8713 

9.8908 

7.7152 

4.8099 

8.3555 
8.789 



4.3149 

4.6747 

4.279 

6.0283 

7.6481 
6.0989 

6.6385 

6.0494 

6.0823 

2.8371 

6.3210 

6.7961 
5.7501 

6.7132 

4.5452 

5.8903 
5.5423 

4.3021 

6.0264 

5.0159 

4.8238 
6.398 



261.69 
207.77 

218.56 

83.213 

140.77 
31.688 

55.119 

93.207 

24.374 

66.192 

70.608 

48.11 
51.422 

90.597 

210.04 

44.337 
55.451 

248.93 

57.745 

148.89 

97.58 
46.158 



0.005 

0.025 

0.015 

0.005 

NA 
NA 

0.015 

0.001 

0.045 

0.015 

0.005 

NA 
0.007 

0.04 

0.015 

0.01 
NA 

0.015 

NA 

0.04 

0.005 
NA 



0.05 
0.140 

0.07 

0.100 

NA 
NA 

0.024 

0.009 

0.01 

0.150 
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NA 
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0.12 

0.05 
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NA 
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NA 
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F2 

NA 

F2 
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NA 

G5 

F8 

F8 

F8 

B9.5V 
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NA 
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201.669 

203.192 
205.697 

211.968 

211.626 

211.207 

212.364 

227.236 

229.203 

232.32 

232.966 

234.948 

234.306 

236.792 

241.475 

241.605 

242.071 

245.118 

252.365 
252.731 

258.477 



17.4714 

16.1046 
-17.759 

18.8854 

11.1521 

15.5078 

18.3606 

15.5816 

11.7905 

13.0392 

14.4351 

15.6862 

12.8062 
-18.335 

10.6215 

16.7817 

18.6567 

17.1348 
-12.911 
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0.2 
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0.25 
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0.2 
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0.2 

0.15 
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0.7 
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2.3347 
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0.883 
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0.3583 
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2.7191 
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2.1077 
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7.8799 
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8.4074 
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8.7012 

9.5831 

7.2047 

7.0607 

9.1309 

9.6686 

9.726 

7.9329 

9.5769 

9.5983 

8.3735 
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6.2926 

6.0757 
6.7612 
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6.2785 

3.4011 
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NA 
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0.025 

NA 
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NA 
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0.04 

NA 
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NA 
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NA 
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NA 
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NA 
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266.238 
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269.851 

273.948 

272.64 

273.655 

275.454 

274.259 

277.105 

276.160 

282.204 

282.714 
283.613 
286.491 
288.756 
294.679 
295.571 
298.277 
297.611 



18.5839 
13.5493 
14.7258 
17.3986 
18.4224 
16.7493 
13.4002 
16.3246 
15.5087 
16.7012 
13.6748 
18.6322 

13.9092 
12.6098 
19.4813 
15.9766 
17.0942 
14.9705 
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